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ABSTRACT 

X-rays produced by compact flares co-rotating with a Keplerian accretion disc 
are modulated in time by Doppler effects. We improve on previous calculations of 
these effects by considering recent models of intrinsic X-ray variability and compute 
the expected strength of the relativistic signal in current data of Seyfert galaxies and 
black hole binaries. Such signals could clearly be seen in, for example, recent XMM- 
Newton data from MCG-6-30-15, if indeed the X-rays were produced by co-rotating 
flares concentrated toward the inner disc edge around an extreme Kerr black hole. 
Lack of the signal in the data collected so far gives support to models, where the 
X-ray sources in active galaxies do not follow Keplerian orbits close to the black hole. 
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1 INTRODUCTION 

Matter accreting onto compact objects (black holes, neu- 
tron stars) can be accelerated to relativistic velocities. The 
best observational evidences of such velocities come from 
distortions of spectral features due to Doppler energy shifts. 
One of the best evidences for relativistic motion in accret- 
ing black holes (both in Active Galactic Nuclei, AGN, and 
black hole binaries, BHB), are the broad Ka lines of Fe near 
6.5 keV (review in Reynolds & Nowak 2003). These can be 
explained by Doppler and gravitational energy shifts, from 
Keplerian motion of the emitting plasma, deep in the grav- 
itational potential wells of the central objects. 

Spectral distortions are not the only effect of relativistic 
motions of sources of radiation. The observed flux should 
be modulated in time, corresponding to periodic Doppler 
enhancement and reduction as the source revolves around 
the central body. The primary emission will be modulated, 
if its source also participates in the rapid Keplerian mo- 
tion. We can thus expect periodic time modulation of the 
primary X-ray emission from accreting compact objects in 
those spectral states, when the emission is coming from mag- 
netic flares. This applies to, for example, BHB in soft states. 
Generally, such time modulation should be observed from 
sources showing broad Fe lines, again, under the assump- 
tion that the primary emission originates in magnetic flares 
rotating with the disc, rather than, e.g., a base of a jet. The 
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presence of a magnetically driven active corona is naturally 
expected from evolution of magnetic field inside an accre- 
tion disc (e.g. Galeev, Rosner & Vaiana 1979), and it is in 
fact supported by recent magnetohydrodynamical simula- 
tions (Turner 2004). 

The purpose of this paper is to estimate the strength of 
the modulation in X-ray power density spectra (PDS) from 
accreting black holes. First computations of the appearance 
of a source orbiting a black hole were performed by Cun- 
ningham & Bardeen (1972) and Bardeen & Cunningham 
(1973) The influence of relativistic effects on time variabil- 
ity of accreting black holes was estimated by Abramowicz 
et al. (1991) and subsequent papers, e.g., Bao (1992), Bao 
& Ostgaard (1995). They assumed a collection of uncorre- 
cted emitting hot spots, so the variability came in part from 
the emergence and disappearance of the hot spots, and in 
part from the relativistic effects. We are more specific in 
our approach as to the form of the intrinsic (i.e. before the 
relativistic effects) variability. 

We assume that the basic form of variability PDS is 
realized without the relativistic effects (we use the flare 
avalanche model of Poutanen & Fabian 1999; see below). 
The latter give additional observable results. The motiva- 
tion for formulating the basic model without relativistic ef- 
fects is that the observed PDS has the basic form of (at least 
approximately) a doubly broken power law for many kinds 
and spectral states of accreting compact objects. This in- 
cludes both the low/hard and high/soft states of black hole 
X-ray binaries, where the accretion flow geometries seems 
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to be quite different, yet the power spectra are similar, at 
least in their high-frequency parts (reviews in Done 2002, 
Zdziarski & Gierlihski 2004). The X-ray power spectra of 
Seyfert galaxies have the same form, with break frequen- 
cies roughly scaling with central mass (Czerny et al. 2001; 
Markowitz et al. 2003). Also, the time scales of broad band 
X-ray variability in accreting compact systems are much 
longer than the Keplerian time scales of the inner accre- 
tion disk, so other mechanisms of generating the variability 
must operate. We note that Miniutti & Fabian (2004) at- 
tribute much of the complex X-ray variability of MCG-6- 
30-15 specifically to Doppler and light bending effects from 
close to the central black hole. At the same time, Vaughan, 
Fabian & Nandra (2003) noticed similarity of power spectra, 
time lags and coherence function between MCG-6-30-15 and 
Cyg X-l, with a simple scaling with the black hole mass. 

However, even though the basic variability model does 
not seem to require relativistic effects, these effects are 
quite clear in energy spectra of many sources. In extreme 
cases, model fits to the observed profiles of Fe Ka line in- 
dicate that the emission must come from well within 6R S , 
the marginally stable orbit in Schwarzschild metric. This is 
claimed to be the case for Seyfert 1 galaxy MCG-6-30-15 
(Fabian et al. 2002), black hole binaries GX 339-4 (Miller et 
al. 2004) and XTE J1650-500 (Miller et al. 2002), where the 
estimated inner radii are « 2-3R g . In addition, the radial 
cmissivity of the line seems to be strongly peaked toward 
the inner radius. Thus, the relativistic effects should con- 
tribute to the variability characteristics, for example, the 
power spectra. 

We adopt the flare avalanche variability model of Pouta- 
nen & Fabian (1999; hereafter PF99), and assign radial lo- 
cations to the flares. We implement the relativistic effects in 
extreme Kerr metric (a = 0.998), and compute the modu- 
lation in the observed X-ray light curves due to relativistic 
effects, as a function of inner disc radius and inclination an- 
gle. We compute the PDS, and demonstrate that for the 
extreme cases of inner radius being close to the marginally 
stable orbit (J? ms = 1.237? g ), the signal should be observed 
in current high quality XMM data. 



2 THE MODEL 

2.1 Variability model and energy spectra 

The model is very similar to that used by Zycki (2002), 
and it is based on the variability model of PF99. The X-ray 
emission is assumed to be produced in active region filled 
in with hot plasma, related to, e.g., magnetic activity in 
accretion disc corona. An active region produces a flare of 
radiation by inverse-Compton upscattering of soft photons 
from the accretion disc. The soft photons, in turn, come 
mostly from reprocessing of the hard X-rays. The spectrum 
of emitted radiation evolves in time during the flare, from 
softer to a rather hard one close to the flare peak. This 
evolution can be caused by, e.g., rising the active region 
above the disc during the flare. The flares are correlated 
to produce longer "avalanches", which help explaining the 
/ _1 part of the power spectrum (see PF99 for a detailed 
description) . Compared to the original formulation of PF99 
we modified the shape of the heating function, lh{t) (where 



I = (L/ Rfiarc){o~T /m c c )) is the compactness parameter). 
It is assumed to be a double exponential function with the 
rise time much longer than the decay time, r r 3> r<j, based 
on the auto- and cross-correlation analysis by Maccarone, 
Coppi & Poutanen (2000). The rise time scales of the flares 
are distributed between r m i n and r max with the distribution 
Pr(r) oc t~ v (see PF99). The cooling of the hot plasma is by 
soft disc photons from reprocessed X-rays, l s (t) = lh(t)D(t), 
and the feedback function, D(t), is defined as in PF99. 

We note that the model of PF99 was originally formu- 
lated for hard/low state of BHB, where the PDS has the 
characteristic shape of a doubly broken power law. Power 
spectra in the soft state have somewhat different form, since 
the low frequency break (from f° to f^ 1 ) is not observed 
(Cyg X-l, Churazov, Gilfanov & Revnivtsev 2001; NGC 
4051 M c Hardy et al. 2004), at least down to frequency 
?s 10~ 5 x /h,br, where / hibr is the high break frequency. Also, 
the shape of the break in the soft state PDS is more grad- 
ual than in hard state. However, these differences do not 
affect our reasoning and computations, because the consid- 
ered effects appear on the f~ 2 part of PDS, i.e. past the 
high-frequency break. This part of PDS is formed by su- 
perposition of Fourier transforms of individual flare profiles, 
and as such, depends on properties of individual flares and 
distribution of their life times. We also note than in the ex- 
treme soft state of Cyg X-l the rollover frequency seems to 
be similar to that in the hard state (fig. 4 in Zdziarski & 
Gierlihski 2004), i.e. there seems to be no difference in the 
PDS extent between the hard and soft states in Cyg X-l. 

Flares are assumed to follow Keplerian orbits, starting 
from a randomly selected azimuthal location, so that the 
primary emission is modulated by the relativistic and kine- 
matic effects. The vertical velocity of a rising plasma region 
(PF99), v z , is usually much smaller than the Keplerian ve- 
locity, 



(where H is the height above the disc of the point of flare's 
maximum luminosity) since H/R<^. 1. Therefore we neglect 
the vertical motion while computing the Doppler effects. For 
the purpose of computing the relativistic effects we assume 
that the emission comes from the equatorial plane. 

Crucial model assumption is the radial distribution of 
flares. A steep centrally peaked distribution obviously en- 
hances the importance of relativistic effects, while a flatter 
distribution makes those effects less important. We assume 
here that the flares are distributed in such a way that the 
energy emitted in X-rays follows the prescription for grav- 
itational energy release in extreme Kerr metric, Q(r) (see, 
e.g., Krolik 1999, p. 152). The steep emissivity is also mo- 
tivated by the results of fitting the observed profile of the 
Fe Ka line, which in some cases indicate very steep and 
centrally peaked emissivity (Fabian et al. 2002; Miller et 
al. 2002, 2004). The above assumption of emissivity is an 
important change compared to our previous model (Zycki 
2002). There the assumption was made that the radial loca- 
tion of a flare is related to its time scale, r tx r'^ z derived 
from th tx , as expected from numerical simulations (e.g. 
Romanova et al. 1998) . This however leads to contradictions 
with observations, since, with all the assumptions of the orig- 
inal model of PF99, the longer flares emit more energy than 



© 0000 RAS, MNRAS 000, 000-000 



Relativistic effects in X-ray variability 3 



the short ones, yet they are located farther away from the 
center (see Zycki 2002 for detailed discussion) . Therefore, in 
the present model there is no relation between the flare time 
scale and its location, but we also discuss in Sec 13. II if relax- 
ing this assumption can modify our results. Radial position 
of a spontaneous flare is generated from probability distribu- 
tion Prif) oc rQ(r) (corresponding to energy dissipation per 
ring), and all the flares stimulated by a given spontaneous 
flares are located at the same radius. 

At each moment of time the energy spectrum is com- 
puted. We approximate the Comptonization spectrum as a 
sum of a power law with exponential cutoff and the Wien 
spectrum, following Zdziarski (1985 and references therein). 
The spectral slope and electron temperature are computed 
using formulae of Beloborodov (1999a, b). The reprocessed 
component is not computed. The Comptonized spectrum, 
A?em, is normalized to lh- 



2.2 Effects of photon propagation in Kerr metric 

We consider a Kerr black hole characterized by its mass, 
M, and angular momentum, J. We use the Boyer-Lindquist 
coordinate system (t, R,8,<f>). The following dimensionless 
parameters are used below 



R_ 

Re 



Ct_ 



J 



dt 



where R g — GM/c 2 is the gravitational radius. We assume 
the highest expected value of the spin parameter, a = 0.998 
(Thorne 1974). 

We find the photon number flux observed from a source 
located at r em , iVobs(rem), by integrating the photon num- 
ber flux in the source rest frame, N eul , over all initial direc- 
tions for which photons reach the observer, and correcting 
the photon emission rate by the time dilation factor, dr/dt, 
where dr and dt are time intervals measured in the rest 
frame and at infinity, respectively (see, e.g., Bardeen 1973). 

All the photon propagation effects in our model are 
treated by the transfer function, T, defined by 



N obs {E,t,i, 



dr 



T(r em ,i,g,A£)N em (E/g,t- At)— dp, (2) 



where i is the inclination angle, At is time interval be- 
tween emission and detection of a given bunch of pho- 
tons, At — AtR g /c, and g is the shift of photon energy, 
9 = E b B /E cul . Note that the transfer function is indepen- 
dent of M. The only dependence on the black hole mass 
occurs in physical scaling of photon travel time and disc 
area. Then, total observed spectrum is given by 



N tot (E,t,i) 



(£)' 



N obB (E,t,i,r)Q(r)2nrdr, 



(3) 



where function Q(r) gives the radial distribution of the emis- 
sion, ri n and r ou t are the inner and outer radius of emitting 
region (r ou t is fixed at 100i? g ) and D is the distance to the 
source. 

We construct the transfer functions following a Monte 
Carlo method described in detail in Gierlihski, Maciolek- 
Niedzwiecki & Ebisawa (2001). However, Gierlihski et al. 



(2001) consider time-averaged emission, while our model 
treats time-dependent spectra. Therefore, we extend the 
procedure to include additional information on the time of 
photon propagation, as follows. 

A large number of photons are emitted from a point 
source on a Keplerian orbit in the equatorial plane. The 
angular velocity of the source with respect to the distant 
observer 



r 3 / 2 + a 



(4) 



We assume that the source is isotropic and then initial direc- 
tions of photons are generated from the distribution uniform 
in both </>em and COS Pern i 

where 6 cm is a polar angle between 
the photon initial direction and the normal to the equato- 
rial plane and om is the azimuthal angle, in the equatorial 
plane, with respect to the r-direction. Solving equations of 
photon motion, as described below, we find the energy shift, 
inclination and arrival time, (g, i, At). 

Each element of the transfer function, T(r, i, g, At), is 
computed by summing all photon trajectories emitted from 
r for all angles (<f> eul , e m) for which required (g, i, At) are 
obtained. 

The travel time and the change of azimuthal angle of 
a photon emitted at a distance r from the equatorial plane 
are given by (e.g. Misner, Thorne & Wheeler 1973) 



, A — a sin 2 8,„ f a(r' 2 + a 2 — Xa) . ._. 
A ^ = / SVW? d * + / A^ 2 (5) 



vr/2 



(r 2 + a 2 )(r 2 + a 2 - Xa) 
ASK 1 / 5 



dr, (6) 



tt/2 



where K and are radial and polar effective potentials, 
respectively, 



&(r) = (r 2 + a 2 - Xa) 2 — A [£ 2 + (A — a) 2 ] 
6(0) = i 2 + cos 2 6 (a 2 - A 2 /sin 2 0) , 



(7) 



£ and A are photon constants of motion, related to the emis- 
sion angles by 



A 



+ V 



cos 6 C 



( r 2 A i/2 + 2arV)/A + Q K sin6U sin 

£ = (^) 1/2 (i-^ 2 r 1/2 (i-Afi K ) 



A = r 2 - 2r + a 2 

A = (r 2 + a 2 ) 2 -a 2 A 

V = {n K -2ar/A)- 



A 



•a A 1 / 2 : 



(8) 

(9) 
(10) 

(11) 



and the angle, i, at which the photon is observed far from 
the source is determined by the integral equation of motion 



5T 1/2 dr= / 6~ 1/2 d<9. 

Jtt/2 

The shift of energy is given by 
A\V2 (i~V 2 ) 1/2 



g = r 



(Ay (i 
\a) ~ 



- n K x 



(12) 



(13) 
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A photon with trajectory characterized by A0 p h must 
be emitted when the source is located at <f> 3 — <f) hs — A(/> P h 
to reach the observer whose azimuthal angle is o bs- Then, 
taking into account c/> s (t) — (j>o + Qxt, where <f>o is a random 
initial location, we find the arrival time of the photon At = 
Af ph - A(/> ph /Q K for A</>p h < and At = Af ph - (2ty - 
A(j> ph )/Sl K for Acj> ph > 0. 

The time dilation factor is given by 

The transfer function for a given radius is normalized to 
the total number of emitted photons. This means that 
in includes the effect of photon capture by the black 
hole/accretion disc. This effect suppresses the emission from 
close to the black hole reaching a distant observer. 

Our transfer function is not ideally smooth when plot- 
ted, for example, as a function of A</> p h for a given radius, 
which is a consequence of finite statistics in the Monte Carlo 
method. This is apparent in some cases, at highest frequen- 
cies of PDS. 

The sequence of spectra created by the above proce- 
dure is subject to standard analysis in the time and Fourier 
domains (see e.g. van der Klis 1995; Nowak et al. 1999; 
Poutanen 2001). The power spectra are computed from light 
curves at 3 keV. 

All computations are done assuming the central black 
hole mass M = 10 6 Mq. The parameters of intrinsic vari- 
ability model are assumed to follow a simple scaling with 
the black hole mass. That is, their values are adopted as for 
Cyg X-l and then scaled by a factor of 10 s 



3 RESULTS 

The observed flux of radiation from a flare co-rotating with 
the disc is periodically modulated. Main contribution to the 
modulations comes from the Doppler boost: the intensity is 
modulated as the observed energy of photons varies with 
their point of emission, so that the invariant Iv/v 3 is con- 
stant. Additional effect comes from the spectral shape: since 



N h s {E) contains contributions from N exn (E/g) oc g F E~ F 
(where V ~ 1.7 is the photon spectral index), this en- 
hances the Doppler modulation. Finally, there is also pe- 
riodic modulation due to varying At p h, which is influencing 
the times/rate of arrival of photons. 

In order to demonstrate the magnitude of the modula- 
tion, we plot in Fig. the light curves and power spectra 
from an intrinsically constant signal emitted at 1.25 R s and 
6 7? g , and observed at i = 60°. The modulation is very strong 
and non-sinusoidal. For 7? e m = 1.25 7? g its amplitude is 
» 2xl0 3 . The non-sinusoidal shape of the modulation means 
that the PDS contains strong harmonics - peaks at multiples 
of the fundamental frequency, /o = = (27r/f2K) -1 - The 
case of -Rom = 6i? g is qualitatively similar, with the funda- 
mental frequency smaller by factor Qk(G -R g )/f2K(1.25 R g ) ~ 
6.5. The underlying PDS continuum is approximately a 
power law P(f) oc f~ 2 , which is a power spectrum of a 
constant signal for t > 0. 

Our main results are presented in Fig. [5] The underly- 
ing PDS breaks from oc / _1 to oc f~ 2 at / br ~ 10~ 5 Hz (for 
10 6 Mq black hole). Superimposed on this continuum is a 
strong signal from relativistic modulation. The signal is far 
from quasi-periodic, both because of the assumed range of 
radii where the flares are located, and because of the strong 
harmonic content. For R ln = 6i? g the orbital frequency at 
the inner radius can actually be seen in the PDS, together 
with two of the higher harmonics. The feature drops sharply 
at / = T~l(R in ) » 2 x 10" 3 Hz, while its slope at low / 
side is related to the radial distribution of the flares (see 
Abramowicz et al. 1991). For _Ri n = 1.25 R s the relativistic 
signal in PDS is quite broad, and the frequency Qk(t = i?i n ) 
is not evident there. The main reason for this is again the 
strong harmonic content, but there are some additional fac- 
tors. Firstly, the energy dissipation rate dependence on ra- 
dius, Q(r), goes to zero at r = J? ms » R[ n . Secondly, both 
the correction for time units, eg. 1141 and the fraction of pho- 
tons escaping to infinity are functions of radius, strongly 
suppressing the emission from small radii. All these com- 
bine to produce the luminosity received from a given radius, 
L(r), which is peaked at somewhat larger radius than Ri n 
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Figure 2. Power density spectra from the flare model, with X-ray flare emission modulated by relativistic effects from Keplerian motion 
of the flares, in the Kerr metric with a = 0.998. The flares are concentrated toward the inner disc radius, so that the radial emissivity 
follows that expected for an accretion disc in the same metric. Inner disc radius is R ln = 1.25 R g (a), and R ln = 6 R g (b). Results for a 
range of inclination angle, i, are plotted. The PDS for i = 90° is combined from two: one from a light curve with tbi n = 10 s, the other 
with tt,in = 1 s. The thin magenta histograms in both panels show the PDS without the relativistic effects. Black hole mass is 10 6 Mq . 
For a BHB system with 10 Mq the additional signal appears above 100 Hz. 



(the peak is at « 2R S ). Also, the radial function L(r) is 
somewhat flatter than r -3 . 

There is a larger spread in the overall PDS normaliza- 
tion for Ri n — 1.25R g than for R ln — 6R g . This is most 
likely caused by smaller fraction of flares (thus larger dis- 
persion) contributing to the observed emission for smaller 
radius. The reason is again that the flares are emitted ac- 
cording to the local energy dissipation rate, Q(r), while in 
the observed emission the contribution from smallest radii 
is relatively suppressed. 

In the considered frequency range, 10 -J < / < 0.05 Hz, 
the integrated r.m.s. is m 0.09 without the relativistic ef- 
fects. The additional integrated signal contribution to r.m.s. 
depends, obviously, on Ri n and i. For Ri n — 1.25 R g it is 
0.03-0.40 (see Fig. 01; specifically, it is « 0.06 for i = 30°. 
For Ri n — 6 R s the additional variability r.m.s. is in the 
range 5 x 10" 3 -0.11. 

For the assumed mass of 10 6 Mq , the additional signal 
appears between 10 -3 and 10 -2 Hz, i.e. on timescales of 
hundreds of seconds. This would correspond to / > 100 Hz 
for a BHB with 10 M . 

3.1 Dependence of the strength of the relativistic 
signal in PDS on radial distribution of flare 
time scales 

The strength of the relativistic signal in PDS may depend 
on assumptions about possible correlations between radial 
position of a flare and its time scale. In our basic model 
there is no such relation, which means that the probabil- 
ity Pt-(th) is independent of radius. The relativistic signal 
is reduced if longer flares occur preferentially at larger dis- 
tance, because - at a fixed modulation period - the shorter 
the modulated flare, the lower the quality factor of result- 
ing periodic modulation. We repeat, however, that the sim- 
ple physically motivated prescription uniquely linking r and 
Tfl, Tfi oc HkM -1 leads to incorrect radial distribution of 



emitted energy, which increases with radius (Zycki 2002). 
We have therefore tested more complex j — r n relations. We 
modified the radial probability distributions, multiplying the 
original function, P r (r), by a gaussian, exp[— (r — ro) 2 /a 2 ], 
where the center ro oc r^f 3 . We have made a number of trial 
runs with various values of a, either constant, or dependent 
on ro. This form of modified probability distribution, P r (r), 
indeed favours longer flares to appear farther away from the 
center, and the strength of this effect can be controlled by 
cr. For a — > do, P r (r) reduces to original rQ(r), while for 
a — > 0, the unique relation ro oc r^f 3 is recovered. Obvi- 
ously, the above form of P r (j) is completely arbitrary; its 
only purpose is to test the dependence of our results on pos- 
sible unknown relation r-r u . 

Fig. shows examples of PDS and radial distributions 
of emission for such modified probability distributions. The 
red histograms show a case of strong correlation between 
r and tr, with a — 0.2ro. This results in energy emission 
shifted towards larger radii, and, in consequence, no rela- 
tivistic signal in PDS. The green histograms show results 
for a/r = 0.2 + 0.4 ln(r fl /r min )/ ln(r max /r min ), meaning a 
narrow radial distribution for short lived flares, which pro- 
gressively becomes broader for longer flares. This gives the 
emissivity somewhat shifted to larger r and weaker QPO 
signal compared to the basic unmodified model. However, 
when the function P r (r) is now modified to be a steeper 
function of r, so that the radial emissivity is restored to its 
original dependence, rQ(r), the strength of the QPO signal 
is also restored (the flares are moved closer to the center). 

Concluding, when the radial emissivity is forced to have 
a prescribed form, the (unknown) relation between radial 
position and life time of a flare does not seem to have any 
influence on the strength of the relativistic signal in PDS. 
We do note that this requires that some longer flares occur 
close to the center, in apparent contradiction with the ba- 
sic condition that the flare time scale should be of order of 
Keplerian time scale (e.g. Romanova et al. 1998). 
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Figure 3. Examples of modified radial probability distribution of flares, with resulting radial emissivities (a) and PDS (b). In the 
modified models longer flares occur preferentially at large distance, as expected from physical arguments (see Sec IS. II , The emissivities 
were multiplied by r 3 (asymptotically, Q(r) oc r~ 3 ). Black histograms show original model for Ri n = 1.25 R g and % = 30° as a reference. 
Red histograms show a model with strong correlation between r and Tfl, resulting in the emissivity peaking at large r and lack of the 
relativistic signal in PDS. Green histograms demonstrate a case of relatively narrow radial distribution of short lived flares, while the 
radial distribution of longer flares becomes progressively broader. When the radial emissivity is forced to have a given steep radial 
dependence, the strength of the relativistic signal in PDS is independent of any r-Tfl relation. 



4 DISCUSSION 

The broad Fe Ka lines and reflection spectra observed in 
X-ray data of many accreting black holes imply fast Kep- 
lerian rotation of the reprocessing matter, but they do not 
imply that the sources of primary radiation participate in 
this motion. Moreover, the primary spectra, being of an ap- 
proximately power law form, are not expected to be dis- 
torted by the relativistic effects. Therefore, we have com- 
puted timing signatures of these relativistic effects on pri- 
mary Comptonized emission. We have assumed that the ba- 
sic form of variability PDS is a result of a flare avalanche 
model of PF99 and that the flares correspond to compact 
sources co-rotating with the Keplerian accretion disc. We 
have estimated the position and strength of the additional 
signal in PDS due to periodic modulation of intensity from 
the Doppler and other effects. We now discuss briefly the ap- 
plicability of our results to current data of Seyfert galaxies 
and black hole binary systems. 



4.1 Seyfert galaxies 

The Seyfert 1 galaxy MCG-6-30-15 is the best example of 
an AGN with a broad Fe Ka line. The line profile in this 
source apparently requires X-ray reprocessing to take place 
very close to the horizon of a rapidly rotating black hole 
(e.g. Fabian et al. 2002 and references therein) . If the sources 
of primary emission participate in the Keplerian emission, 
this would produce very clear timing signatures. In order to 
see whether they would be observable in current data, we 
simulate a light curve corresponding to the recent XMM- 
Newton observations described by Vaughan et al. (2003). 
We use the unmodified model, i.e. without any r-rfl corre- 
lation (Sec. l3.l1 . since it ensures the steep radial emissivity, 
as required by the Fe line profile fits. The model light curve, 
F(t), is scaled to a given mean value, fi. "Observed" light 



curve is then formed by drawing a number of counts in each 
time bin from Poisson distribution of mean F(t n ), where 
t D = n X thin, n = 1,2, . . . , N . Thus simulated light curve 
contains the effect of Poisson noise. The time bin is tbin = 10 
s, total duration of the light curve is » 10 5 s, while the mean 
value, n = 20 cts/s, all as in actual observations (see fig. 1 in 
Vaughan et al. 2003). Fig. 0] shows the resulting PDS with 
the Poisson noise level subtracted, together with a best fit 
power law plus Lorentzian model. A simple power law model 
gives a bad fit to the data, xt = 747/61 d.o.f. Improvement 
by adding the Lorentzian is significant, = 257/58 d.o.f., 
although the description of the additional feature in PDS is 
not unique: a broken power law model gives a fit of the same 
quality. We have checked that similarly significant difference 
between the two model fits to the PDS with relativistic ef- 
fects is maintained, even if the the mean source count rate is 
4 time lower. Thus, the signal from relativistic modulation, 
if present, could be observable, in current data. The signal 
is, however, not observed (see fig. 4 and fig. 5 in Vaughan et 
al. 2003) . The lack of the signal is problematic for the model 
of rotating flares, suggesting alternative geometries, which 
do not require Keplerian motion of the X-ray sources. 

On the other hand, periodically variable redshifted Fe 
Ka line was recently reported in XMM-Newton data of NGC 
3516 by Iwasawa, Miniutti & Fabian (2004). The variabil- 
ity could be interpreted as a result of Doppler effect from a 
flare co-rotating with a disc, with orbital period ~ 25 ksec. 
Since this kind of transient modulated line emission is not 
expected from a model with axially symmetric illumination, 
this observation provides an argument in favour of the ro- 
tating source model. However, the primary continuum does 
not show the expected modulation. These two observational 
results - modulation of Fe Ka line and lack of modulation 
of primary continuum - if firmly established, will certainly 
require revision of current X-ray emission models. 

One obvious consequence of the relativistic modulation 
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Figure 4. The relativistic signal could be easily seen in PDS from high quality XMM-Newton data of MCG-6-30-15. The model light 
curve for _R; n = 1.25 R g and i = 30° (see Fig. was used to generate a light curve of quality corresponding to the observations of 
Vaughan et al. (2003). Total length of the light curve is fa 4 X 10 5 s, the mean count rate is 20 cts/s, and the time bin is 10 sec. Panel (a) 
shows the resulting PDS with Poisson level subtracted, and with a best fit power law plus a Lorentzian model. The right panels show the 
data— model residuals. The additional signal in PDS is highly significant, although its detailed description is not unique. No such signal 
is seen in the actual data of that Seyfert galaxy (see fig. 4 and fig. 5 in Vaughan et al. 2003). 



is that highly inclined sources should show stronger modu- 
lation than less inclined ones. This was also pointed out by 
Czerny et al. (2004) for Seyfert 1 vs Seyfert 2 galaxies. This 
effect may not be easily observable, however, if the geometry 
of accretion in Seyfert galaxies corresponds to low state of 
BHB, and the accretion discs are truncated far from the last 
stable orbit. Since the strength of the relativistic signal de- 
creases very quickly with increasing inner disc radius, these 
effects can only be observed in sources in high states where 
the disc extends all the way to the last stable orbit. 



5 CONCLUSION 

Keplerian motion of sources of X-ray radiation in accreting 
black holes gives a very strong signature in power spectra, 
if the sources are located close to the center. Such signal, if 
present, would be observed in recent high quality data, for 
example the XMM-Newton data from the Seyfert 1 galaxy 
MCG-6-30-15. Lack of the signal in those data posses prob- 
lems for models assuming that the X-ray sources participate 
in the Keplerian disc motion, as for example expected for 
magnetic flares. 



4.2 Black hole binaries/other objects 

Two BHB were reported to have very broad Fe Kcv lines in 
their XMM-Newton spectra: GX 339-4 (Miller et al. 2004) 
and GRO J 1650-500 (Miller et al. 2002). The inner disc ra- 
dius inferred from the fits to the line profiles are 2-3_R g , 
similar to that in MCG-6-30-15. Thus, in principle, high- 
frequency PDS should contain the signal due to relativistic 
modulation. However, as is well known, stellar-mass black 
hole systems are less suitable to investigating short time- 
scale variability than AGN: even though their X-ray fluxes 
are ~ 10 3 times higher than from AGN, the required time- 
scales are 10 5 -10 7 times shorter. Our simulations indicate 
that it would be rather difficult to detect the additional sig- 
nal above 100 Hz. On the other hand, high-frequency quasi- 
periodic oscillations (QPO) are observed in power spectra of 
X-ray binaries. Some of those QPO are located at high fre- 
quencies, / > 100 Hz, i.e. in a region where the signal from 
relativistic modulations might appear (see, e.g., Schnittman 
& Bertschinger 2004). It is rather uncertain, however, how 
this mechanism might explain the observed QPO, since it 
would require modulation from a narrow ring of the disc, 
and suppression of the harmonics, in order to produce a 
narrow QPO feature. 
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